function A = skew(A)
% A = skew(A)
% Returns the skew-symmetric part of A

  A = 0.5 * (A - A');

end